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Abstract. The masses of compact objects like white dwarfs, neutron stars 
and black holes are fundamental to astrophysics, but very difficult to measure. 
We present the results of an analysis of subluminous B (sdB) stars in close 
binary systems with unseen compact companions to derive their masses and 
clarify their nature. Radial velocity curves were obtained from time resolved 
spectroscopy. The atmospheric parameters were determined in a quantitative 
spectral analysis. Based on high resolution spectra we were able to measure the 
projected rotational velocity of the stars with high accuracy. In the distribution 
of projected rotational velocities signs of tidal locking with the companions are 
visible. By detecting ellipsoidal variations in the lightcurve of an sdB binary 
we were able to show that subdwarf binaries with orbital periods up to 0.6 d 
are most likely synchronized. In this case, the inclination angles and companion 
masses of the binaries can be tightly constrained. Five invisible companions have 
masses that are compatible with that of normal white dwarfs or late type main 
sequence stars. However, four sdBs have compact companions massive enough 
to be heavy white dwarfs (> IM©), neutron stars or even black holes. Such 
a high fraction of massive compact companions is not expected from current 
models of binary evolution. 



1. Introduction 

The mass of a star is its most fundamental property. However, a direct measure- 
ment is possible in a small number of binary stars only. White dwarfs, neutron 
stars and stellar mass black holes are the fibal products of stellar evolution. In 
most binaries such faint, compact objects are outshone by their bright compan- 
ions and therefore their orbital motion cannot be measured. Only lower limits 
to the companion mass can be derived. With the analysis method shown here, 
these limitations can be partly overcome. 

Subluminous B stars are considered to be helium core burning stars with 
very thin hydrogen envelopes and masses around 0.5 Mq. Different formation 
channels have been discussed. As it tur ned out, a large frac tion of the sdB stars 
are members of short period binaries ( Maxted et al.l 120011 ). For the se systems 



common-envelope ejection is the most probable formation channel (|Han et al 



20021 ). In this scenario two main-sequence stars of different masses evolve in a 



1 



2 



Geier et al. 



binary system. The heavier one will first reach the red giant phase and fill its 
Roche lobe. If the mass transfer to the companion is dynamically unstable, a 
common envelope is formed. Due to friction the two stellar cores loose orbital 
energy, which is deposited within the envelope and leads to a shortening of the 
binary period. Eventually the common envelope is ejected and a close binary 
system is formed, which contains a helium core burning sdB and a main sequence 
companion. When the latter reaches the red-giant branch, another common 
envelope phase is possible and can lead to a close binary with a white-dwarf 
companion and an sdB star. All known companions of sdB stars in such systems 
are white dwarfs or late-type main-sequence stars. If massive stars are involved, 
the primary may evolve into a neutron star (NS) or a black hole (BH) rather 
than a white dwarf. Since massive stars are very rare, only few sdB+NS or 
sdB+BH systems are expected to be found. 

Since the spectra of most binary sdBs are single-lined, they reveal no infor- 
mation about the orbital motion of their companions, and thus only their mass 
functions can be calculated. 

M^^^pSin^i ^PK^ 
(Mcomp + MsdB)2 27rG • ^ ^ 

Although the RV semi-amplitude K and the period P are determined by 
the RV curve, two of MgdB, ^mmp and sin z remain free parameters. Binary 
population synthesis models fill an et all |2002| ) indicate a possible mass range 



of MsdB = 0.30—0.48 M0 for sdBs in binaries, which underwent the common 
envelope ejection channel. The mass distribution shows a sharp peak at about 
0.46 M0. This theoretical value can be backed up by observations (Edelmann 
these proceedings; Vuckovic et al. these proceedings) as well as asteroseismology 
(e.g. Charpinet et al. these proceedings). 

For close binary systems, the components' stellar rotational velocities are 
considered to be tidally locked to their orbital motions, which means that the 
orbital period of the system equals the rotational period of the companions. If 
the companions are synchronized in this way the rotational velocity v^ot can be 
calculated. 

Vrot = p • (2) 

The stellar radius R is given by the mass-radius-relation. 



R=^^ (3) 

The measurement of the projected rotational velocities Urot sin i therefore 
allows us to constrain the systems' inclination angles i. With Msjb as free 
parameter the mass function can be solved and the inclination angle as well as 
the companion mass can be derived. Because sini < 1, a lower limit for the sdB 
mass is given. To constrain the system parameters in this way it is necessary to 
measure K, P, log 5 and Vrot sinz with high accuracy. 
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2. Observations and Radial Velocity Curves 

Ten stars were observed at least twice with the high-resolution spectrograph 
UVES at the ESO VLT. Additional observations were made at the ESO NTT 
(equipped with EMMI), the Calar Alto Observatory 3.5-m telescope (TWIN) 
and the 4-m WHT (ISIS) on La Palma. Two of the stars (PC 1232-136, 
TONS 183) were observed with the high resolution FEROS instrument at the 
2.2 m ESO telescope at La Silla. The radial velocities (RV) were measured by 
fitting a set of mathematical functions (Gaussians, Lorentzians and polynomi- 
als) to the hydrogen Balmer lines. Sine curves were fitted to the RV data points 
using a minimising method (singular-value decomposition) and the temporal 
power spectrum was generated to obtain a best-fit period (Table 1). 



3. Atmospheric parameters and projected rotational velocities 

The spectra were corrected for the measured RV and co-added. Atmospheric 
parameters were determined by fitting simultaneously each observed hydrogen 
and helium line with a grid of metal-line blanketed LTE model spectra. Partial 
results are listed in Table 1. 



Table 1. Surface gravities, orbital periods, radial velocity semi-amplitudes 
and projected rotational velocities of the visible components. The typical 
error margin for logg is 0.05 dex. f Orbital parameters of this system taken 
from Geier et al. (these pro c eedin gs). J Orbital parameters of this system 
taken from iNapiwotzki et all (|2001f ). 



System 




logs 


P 


K 


Vrot sin i 






[cm^s~^] 


[d] 


[kms^-^] 


[kms^-^] 


HE 0532- 


-4503 


5.32 


0.26560 ± 0.00010 


101.5 ± 0.2 


11.1 ± 0.6 


PG1232- 


-136 


5.62 


0.36300 ± 0.00030 


129.6 ± 0.04 


6.2 ± 0.8 


WD0107 


-3421 


5.32 


0.37500 ± 0.05000 


127.0 ± 2.0 


20.4 ± 0.9 


HE 0929- 


-0424 


5.71 


0.44000 ± 0.00020 


114.3 ± 1.4 


7.1 ± 1.0 


HE0230- 


-4323 


5.60 


0.44300 ± 0.00050 


64.1 ± 1.5 


12.7 ± 0.7 


TONS 183 


5.20 


0.82770 ± 0.00020 


84.8 ± 1.0 


6.7 ± 0.7 


HE 2135- 


-3749 


5.84 


0.92400 ± 0.00030 


90.5 ± 0.6 


6.9 ± 0.5 


HE 1421- 


-1206 


5.55 


1.18800 ± 0.00100 


55.5 ± 2.0 


6.7 ± 1.1 


HE 1047- 


-0436t 


5.66 


1.21325 ± 0.00001 


94.0 ± 3.0 


6.2 ± 0.6 


HE2150- 


-0238 


5.83 


1.32090 ± 0.00500 


96.3 ± 1.4 


8.3 ± 1.3 


HE 1448- 


-0510 


5.59 


7.15880 ± 0.01300 


53.7 ± 1.1 


6.7 ± 2.5 


WD 0048 


-202 


5.50 


7.44360 ± 0.01500 


47.9 ± 0.4 


7.2 ± 1.3 



In order to derive frot sin i, we compared the observed spectra with ro- 
tationally broadened, s ynthetic line profiles. The latter were computed using 
the LINFOR program (|Lemkd[T997l l. Since sharp metal lines are much more 
sensitive to rotational broadening than Balmer or helium lines, all visible metal 
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Figure 1. Companion mass as a function of primary (sdB) mass of the bi- 
naries HE 0532-4503 (left hand panel) and PG 1232-136 (right hand panel). 
The horizontal line marks the Chandrasekhar limit. The dotted vertical lines 
mark the theoretical sdB mass range for the common envelope ejection chan- 
nel (,Han et al.,,20021 



lines were included. A simultaneous fit of elemental abundance and projected 
rotational veloci ty was performed separ ately for every identified line using the 
FITSB2 routine dNapiwotzki et al.ll2004l ). The mean value and the standard de- 
viation were calculated from all measurements. Seeing induced variations in the 
instrumental profile and the noise level were the dominant error sources. Infor- 
mation on the actual seeing conditions for every exposure have been extracted 
from the ESO seeing monitor archive. All other possible sources of systematic 
errors turned out to be negligible. 



4. Nature of the Unseen Companions 

Knowing P, K, log g and f rot sin i we can calculate the mass of the unseen 
companion from Eqn s. 1—3 for any giv en primary mass (see Fig. [T]) . We adopt 
the mass range from iHan et all |2002') , marked by the dotted lines in Fig. [TJ 
There are no spectral signatures of companions visible. Main sequence stars 
with masses higher than 0.45 Mq could therefore be excluded because they would 
contribute to the total flux and could therefore be identified from the spectral 
energy distribution and/or indicative spectral features. The possible companion 
masses can be seen in Table 2. Four of the analysed systems have companion 
masses, which are compatible with either typical white dwarfs (WD) or late 
main-sequence stars (late MS). The companion of WD 0107— 0342 is a massive 
white dwarf. Since the total mass of the binary may exceed the Chandrasekhar 
limit, the system is a good candidate for a double degenerate SN l a progenitor 



with s ub dwarf primary, only the second one after KPD 1930-1-2752 (jGeier et al 
2007al ). 



The very similar HE 0929-0424 and TONS 183 have to have quite mas- 
sive companions. Even in the less likely case that their sdB primaries are of 
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Figure 2. Left panel: The measured Wrotsini of 49 single sdBs is plotted 
against relative fraction of stars as shaded histogram. The solid blank his- 
togram marks the expected uniform distribution of f rot sin i under the as- 
sumption of randomly oriented polar axes and uniform rotational velocity 
Vmt = 8.3 kms^^ for all stars. The dashed blank histogram marks a uniform 
distribution for Uiot = 9.0kms~^ as comparison. All sdBs with iij-ot sinz below 
the detection limit (vertical line) are stacked into the first bin (dotted shaded 
histogram) . The width of the bins is given by the average v^ot sin i uncertainty 
of the sample. Right panel: The measured fj-ot sini of 19 RV variable sdBs is 
plotted against relative fraction of stars as shaded histogram. The solid blank 
histogram marks the uniform distribution derived from the single star sample. 
The width of the bins is again given by the average Wrot sin i uncertainty of 
the sample, which is lower for the RV variable sample due to better quality 
of the data. 

low mass (~ 0.3 M©) the companions would be heavy white dwarfs. At the 
most probable sdB mass, however, their companions would exceed the Chan- 
drasekhar mass limit. There are only two kinds of objects known with such 
high masses and such low luminosities - neutron stars and stellar mass black 
holes. The two systems HE 0532— 4503 and PG 1232—136 have even higher com- 
panion masses, which would exceed the Chandrasekhar limit for any mass of a 
core helium-burning subdwarf star (see Fig. [1]). The three long-period binaries 
HE 1448-0510, HE 2150-0238 and WD 0048-202 could not be solved with the 
described method. These systems cannot be synchronized (see Sect. 5). 

Binaries hosting a neutron star or a black hole are a very rare class of 
objects. From about 50 analysed sdB binaries in our samples, we found two, 
possibly four candidate systems. This fraction of 4-8% is much too high to be 
compatible with any binary evolution model known so far (Podsiadlowski priv. 
comm.). 

Our results depend strongly on the projected rotational velocity. The unex- 
pectedly high masses could only be reduced if the Vmt sin i were underestimated. 
As described above we quantified all possible systematic effects and the overall 
results were very consistent. But even if there would be unaccounted systematic 
effects (e.g. short period pulsations), they would always cause extra broaden- 
ing of the lines. The measured broadening is then due to rotation plus the 
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unaccounted effects, wfiicli means tlie deduced rotational broadening would be 
overestimated. Potential systematic effects yet unaccounted for would therefore 
lead to even higher companion masses. 



5. Orbital Synchronization of sdB Binaries 



Our method rests on the assumption of orbital synchronization. In Fig. [2] we 
compare the f rot sin i-distribution of close sdB binaries to the distribution of 
the single stars. A large fraction of binary sdBs exceeds the maximum Urot = 
8.3kms~^, derived for the single sdB stars, significantly. The most likely reason 
for this is tidal interaction with the companion. 

The question of which mechanisms are responsible for the synchronization 
of such stars is not yet settled. Subluminous B stars have convective cores and 
radiative envelopes. The physical mechanisms fo r tidal dissipation in such stars 
are u nder debate. Different theoretical concepts ( Zahn|[l977l : iTassoul &i Tassoul 
19921 ) predict that the synchronization timescale tsync depends strongly on the 
orbital period but differs in absolute values. Tidal locking can only be assumed if 
the synchronization time scales are significantly lower than the ave rage lifetime 



on the EHB of tEHB ~ 10^ yr. Using the formalism of IZahnI ( 19771 ) igync would 



exceed tEHB for periods longer than -P^ahn ~ 0.4 d, whereas this would be the 
tSsouI ~ 2.0 d if we apply the theory of lTassoul fc Tassoull (fl992). 
But as long as the question of tidal synchroniz ation is not settled, all 



case at P}'"^ 



timescales have to be taken with caution ( Zahn 20051 ). Detailed calculations, 
which take the internal structure of sdBs into account, are not available and ur- 
gently needed. Observational constraints are necessary to guide the development 
of more sophisticated models. 

Up to now such observational constraints are rare. Two short period (~ 
0.1 d) sdB binaries with whit e-dwarf companions s how light variation s due to 
their ellipsoidal deformation ( Orosz Wade 19991: Geier et al.ll2007al ) and are 
therefore most likely synchronized. Randall et al.l ( 20051 ) reported the detection 
of shallow ellip soidal variat i ons in the longer period (~ 0.6 d) sdB-|-WD binary 
PG 0101-^039. iGeier et al.l (|2008l ) verified that the light variations are due to 
ellipsoidal deformation and that tidal synchronization is very likely established 
for PG 0101-1-039 (see Fig. [3|). We conclude that this assumption should hold 
for all sdB binaries with orbital periods of less than half a day. Last but not 
least, van Grootel et al. (these proceedings) carried out an asteroseismological 
analysis of the pulsating sdB in the binary Feige 48 (sdB-|-WD, P = 0.376 d) 
and found it to be synchronized. 

In Table 2 we show which binaries fulfill the empirical or theoretical criteria 
for synchronization. It has to be pointed out that three of our candidate systems 
with massive compact companions are consistent with all of them. 



6. Conclusion 

Out of 12 analysed sdB binaries, five have companion masses compatible with 
white dwarfs of typical mass or late-type main-sequence stars. The properties 
of these systems are in full agreement with binary population-synthesis simula- 
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Figure 3. Ligh tcurve of PG 0101+039 with superimposed models 
(|Geier et al.ll2008f l. The two sohd curves confine the best fit models with 
parameters and associated uncertainties derived from spectroscopy under the 
assumption of orbital synchronization and Mg^B = 0.3 Mq. The two dotted 
curves confine models for MsdB = 0.7 Mq. 

tions. Four binaries have surprisingly high companion masses, which leads to the 
conclusion, that the companions have to be white dwarfs of unusually high mass 
or even neutron stars or black holes. This high fraction cannot be explained 
with current evolutionary calculations. As our analysis assumes synchroniza- 
tion, a better understanding of this process for sdB stars is urgently needed. A 
larger sample of sdB binaries has to be studied with our method to search for 
systematic trends and improve statistics. 

The presence of such a high fraction of heavy binaries in our samples raises 
several questions. Is the formation and evolution of sdB stars linked to that of 
heavy compact objects like neutron stars or black holes? Can sdB stars be used 
as tracers to find more of these exotic objects? Is there a hidden population of 
these objects present in our galaxy? 
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Table 2. Inclination angles, rotational velocities, companion masses and 
possible nature of the unseen companions. The lower companion mass corre- 
sponds to an sdB of 0.3 M0, the higher limit to an sdB of 0.48 Mq. Binaries 
up to orbital periods of 1.3 d can be solved consistently under the assumption 
of synchronization. Empirical and theoretical indications for synchronization 
are given for individual stars. 



System 


i 

[dej 


l] 


^rot 

[kms""'] 


-^comp 

[Mo] 


Companion 


HE0532-4503i'2'3'4 


13 - 


17 


47 


1.40 - 3.60 


NS/BH 


PG1232-136^'2'3'4 


14 - 


19 


25 


2.00 - 7.00 


NS/BH 


WD0107-342i'2'3.4 


37 - 


50 


33 


0.48 - 0.87 


WD 


HE0929-04242'4 


23 - 


29 


18 


0.60 - 2.40 


WD/NS/BH 


HE0230-43232'4 


38 - 


50 


21 


0.18 - 0.35 


WD/late MS 


TONS 183^ 


22 - 


29 


18 


0.60 - 2.40 


WD/NS/BH 


HE 2135-3749^ 


66 - 


90 


8 


0.35 - 0.45 


WD/late MS 


HE 1421-1206'' 


56 - 


90 


8 


0.15 - 0.30 


WD/late MS 


HE 1047-0436^ 


62 - 


90 


7 


0.35 - 0.60 


WD/late MS 


HE 2150-0238'' 










no solution 


HE 1448-0510 










no solution 


WD 0048-202 










no solution 



^Empirical indications for synchronization could be found in subdwarf binaries with 
longer periods by asteroseismological modelling. 

^Empirical indications for synchronization could be found in subdwarf binaries with 
longer periods by detection of ellipsoidal deformation. 

•^Theoretical synchroniz ation timesca les are shorter than the average EHB lifetime ac- 
cording to the theory of lZahnI (|1977D . 

''Theoretical synchronization timescales are shorter than the average EHB lifetime ac- 
cording to the theory of iTassoul fc Tassoull ()1992l ). 
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